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One of the fundamental problems of modern cosmology is to explain the origin of all the matter 
and radiation in the Universe today. The inflationary model predicts that the oscillations of 
the scalar field at the end of inflation will convert the coherent energy density of the inflaton 
into a large number of particles, responsible for the present entropy of the Universe. The 
transition from the inflationary era to the radiation era was originally called reheating, and 
we now understand that it may consist of three different stages: preheating, in which the 
homogeneous inflaton field decays coherently into bosonic waves (scalars and/or vectors) with 
large occupation numbers; backreaction and rescattering, in which different energy bands 
get mixed; and finally decoherence and thermalization, in which those waves break up into 
particles that thermalize and acquire a black body spectrum at a certain temperature. These 
three stages are non-perturbative, non-linear and out of equilibrium, and we are just beginning 
to understand them. In this talk I will concentrate on the preheating part, putting emphasis 
on the differences between preheating in chaotic and in hybrid inflation. 



1 Introduction 

At the end of inflation all the energy density is in the homogeneous zero mode of the inflaton 
field. The Universe is in a vacuum-like state with zero temperature and vanishing particle and 
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entropy densities. The problem of reheating is how to convert al 
a state of thermalized relativistic particles. The original analysis H 
decay of the inflaton into bosons and fermions, as if the inflaton were already an ensemble of 
decoherent particles. Reheating ended when the total decay rate of the inflaton was of the order 
of the expansion rate of the Universe, r ~ H, while the total energy of the inflaton field decayed 
exponentially fast into other particles. As a consequence, the final reheating temperature only 
depended upon T. We understand today that, for certain parameter ranges, there is a new 
decay channel that is non-perturbativep due to the coherent oscillations of the inflaton field, 



which induces stimulated emission of bosonic|| particles into energy bands with large occupation 
numbers. The modes in these bands can be understood as Bose condensates, and they behave 
like classical waves. The backreaction of these modes on the homogeneous inflaton field and the 
rescattering among themselves produce a state that is far from thermal equilibrium and may 
induce very interesting phenomena, such as non-thermal phase transitions S with production 
of topological defects, a stochastic background of gravitational waves P production of heavy 
particles in a state far from equilibrium, which may help GUT baryogenesisH or constitute today 
the dark matter in our UniverseQ These classical waves eventually reach a state of turbulence 
where, hopefully, decoherence will occur and thermalization will follow, although these stages 
are not yet fully understood, either analytically or numerically. 

The period in which particles are produced via parametric resonance is called preheating^ 
The idea is relatively simple, the oscillations of the inflaton field induce mixing of positive and 
negative frequencies in the quantum state of the field it couples to. In the language of quantum 
fields in curved space, creation and annihilation operators mix via Bogoliubov transformations]] 
Ofc = c^kO-k + (3*<^-ki an d with every oscillation of the inflaton field, new particles are produced, 
n/j = (0|aj.afc|0) = \(3k\ 2 - I n the case of chaotic inflation, with a massive inflaton <f> coupled 
to a massless scalar field x, the evolution equation for the Fourier modes, Xk + w^X^ = 0, 
with Xk = a 3 / 2 (t)xk an d w\ = k 2 /a 2 (t) + g 2 4> 2 (t), can be cast in the form of a Mathieu 
equation, with coefficients A = k 2 /4a 2 m 2 + 2q and q = g 2 & 2 /4m 2 , where $ is the amplitude 
and m is the frequency of inflaton oscillations, (j>(t) = <I>(i) sin mi. For certain values of the 
parameters (A, q) there are exact solutions that grow exponentially with time, and each mode k 
belongs to an instability band of the Mathieu equation]] These instabilities can be interpreted 
as coherent "particle" production with large occupation numbers. One way of understanding 
this phenomenon is to consider the energy of these modes as that of a harmonic oscillator, 
E)~ = \Xf : \ 2 /2 + u> 2 \X)~\ 2 /2 = huj-fok + 1/2). The occupation number of level k can grow 
exponentially fast, ~ exp(2/ifcmt) 3> 1, and these modes soon behave like classical waves. 
It is analogous to the well-known mechanism of generation of density perturbations during 
inflation^ The parameter q during preheating determines the strength of the resonance. It is 
possible that the model parameters are such that parametric resonance does not occur, and 
then the usual perturbative approach would follow, with decay rate T. In fact, as the Universe 
expands, the growth of the scale factor and the decrease of the amplitude of inflaton oscillations 
shifts the values of (A, q) along the stability/instability chart of the Mathieu equation, going 
from broad resonance, for q 3> 1, to narrow resonance, g< 1, and finally to the perturbative 
decay of the inflaton. Parametric resonance will stop whenever the inflaton decay is dominated 
by the perturbative decay, qm < T, or when the instability modes are redshifted away from the 
(last) narrow resonance band, q 2 m < H. 



2 Preheating in hybrid inflation 

Until recently, preheating had been studied in chaotic and new inflation onlyP where the end of 
inflation occurs when the rate of expansion is of order the mass of the field, m ~ H, and the 
inflaton starts to oscillate around the minimum of its potential, inducing particle production.]] 
In a recent paperp^l we have studied the case of hybrid inflation, where the end of inflation 
is triggered by the symmetry breaking of another scalar field coupled to the inflaton, and not 



"Fermions can also be parametrically amplified, but their occupation numbers are constrained, rik < 1, by 
Pauli's exclusion principle. m 
'For an introduction to particle production in strong external fields, see Grib et alu m 

c For a recent comprehensive review on preheating after chaotic inflation see KofmanU and references therein. 
d Although production of classical "waves" is more appropriate, we will nevertheless use the terminology of 
"particle" production with large occupation numbers. 



by slow-rollS This fact alone gives hybrid inflation several advantages with respect to chaotic 
inflation. First, we can consider models of inflation at low-energy scales, say the electroweak 
scale, which nevertheless give the correct amplitude of temperature fluctuations in the CMB. 
Second, we may have an effective frequency of oscillation of one or both of the scalar fields being 
much larger than the rate of expansion. As we will see, this induces a very emcient and long- 
lived narrow resonance, where many oscillations occur in one Hubble time, and a large number 
of particles are produced before their momenta are redshifted by the expansion of the Universe. 
Preheating can be very emcient in this case. Third, at the stage of rescattering, gravitational 
waves are produced with wavelengths at most of order the size of the horizon at that time, large 
enough today to be detected at gravitational wave interferometers such as LIGO. 

There are two fields in hybrid models, <j) drives slow-roll inflation and a triggers its end 
through a symmetry-breaking potential, 

a) = ^(M 2 - Aa 2 ) 2 + \m 2 <? + ^ W . (1) 

These fields may couple to yet another one, x, with couplings (h 2 <fi 2 + h?, <r 2 )x 2 /2. Preheating 
depends very strongly on the model parameters. During inflation, the symmetry-breaking field a 
has a large mass, due to its coupling to <p, and is fixed at a = 0. The inflaton eft slow-rolls down its 
effective potential V(<p) = Vq + m 2 (j) 2 /2, driving inflation and producing the metric fluctuations 
that later will give rise to temperature anisotropies in the CMB and density perturbations for 
large-scale structure. The amplitude of temperature anisotropies on the scale of the horizon, as 
seen by COBEpj constrains the parameters of the model to satisfy 

-V 9^ - 3.5 xl(T 5 , 2 
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A m 2 Ml 

< 0.2. (3) 

Contrary to the case of chaotic inflation, these constraints leave plenty of freedom to chose the 
model parameters. One of the advantages of hybrid inflation is that the rate of expansion at the 
end of inflation, H ~ M 2 /\/AMp, could be in a wide range of scales, from just below Planck 
scale, as in models of supergravity hybrid inflationjlJ all the way to the electroweak scale. In 
terms of the rate of expansion, COBE constraints can be written as g = 2 x 10 -4 (n — 1) M/H 
and n - 1 = 2m 2 /3H 2 < 0.2. 

When the inflaton <j) falls below (p c = M/g, the field a triggers the end of inflation via 
a sudden or "waterfall" spontaneous symmetry breaking. This process occurs because of the 
exponential growth of quantum fluctuations. In most models of hybrid inflation, the mass- 
squared of the a field changes from large and positive to large and negative in much less than 
one e-fold|] It can then be shown that all tachyonic modes grow at approximately the same speed 
and thus the a field can be described as a homogeneous fieldS The behaviour of (j) and a after 
the end of inflation therefore follows the homogeneous field equations, and depends crucially on 
the ratio of couplings, 

^(n-lAo- 4 . (4) 

Depending on whether g <C A, g 2 3> A or g 2 ~ A, the oscillations around the minimum will 
occur mainly along the (j> field, the a field, or both, respectively. Due to their mutual couplings, 
the resonant .production of particles could be very efficient or completely suppressed, depending 
on this ratioll3 Another important factor is the relation of the frequency of oscillations to the 
rate of expansion. In chaotic inflation this was fixed by the end of slow-roll condition, m ~ H. 

e There are hybrid models with two stages of inflation where this condition is not satisfied 



However, in hybrid inflation we have the freedom to lower the rate of expansion with respect 
to the mass M and thus we can consider models where the field oscillates many times, as much 
as 10 8 times, in one Hubble timefij Thus, the resonant modes remain for a long time in their 
instability bands and may produce large occupation numbers before they are redshifted away 
by the expansions In particular, in hybrid preheating, the narrow resonance may be long-lived, 
q 2 fh > H, even for very small parameters g C 1, thanks to ra » ff. This is a fundamental 
difference with respect to chaotic inflation. Let us now consider in some detail the three different 
cases mentioned above. 

In the first case, g 2 <C A, the frequency of oscillations of the (j> field is fh ~ gM/yX 3> m 
and, depending on the model parameters, can be much greater than the rate of expansion. 
For fh ~ H, the symmetry-breaking field soon settles at its minimum <7o, giving large masses 
to any field it couples to, while the amplitude of inflaton oscillations quickly decays with the 
expansion of the Universe. As a consequence, there is essentially no production of either <p or 
a particlesuJ However, if a massless scalar (or vector) field x couples to 0, and not to a (for 
whatever symmetry reason), it is possible to produce x-particles in the broad resonance regime, 
for a wide range of couplings. This is the situation that most resembles that of chaotic inflation; 
the growth parameter is \i ~ 0.13 and backreaction sets in after about 20 oscillations of the <j> 
field. In the other case, fh ^> H, the amplitude of oscillations of 4> after the end of inflation 
remains large, of order <f) c , during many oscillations, driving explosive particle production even 
in the narrow resonance, for fi^ <C 1. Since the rate of expansion is so small, the modes remain in 
the narrow resonance for a long time, and backreaction on the inflaton oscillations occurs before 
the modes are redshifted away from the resonancet3 This is very different from the behaviour 
in chaotic inflation. However, even if \ particles are quickly produced, we still do not have a 
significant production of <f> or a particles before x-P r °duction backreacts on <f>. 

In the case g 2 3> A, the situation described above is reversed: it is now the symmetry- 
breaking field that oscillates around its minimum, while the inflaton field becomes negligible. 
The frequency of a oscillations is fh = \[2M and, as before, it can be much larger than the rate 
of expansion. For fh ~ H, the inflaton settles at <p = while a oscillates with ever-decreasing 
amplitude. In this case, there is an insignificant production of either cj> or a particles, due to 
their large effective massesH Furthermore, if a field \ couples to either or a its production 
via parametric resonance is suppressed, because of the small amplitude of cf> oscillations or of 
the effective mass induced by uq, respectively. It follows that preheating in this case is very 
inefficient. On the other hand, for fh ^> H, explosive particle production occurs for all fields, 
but mainly for a particles. Even though the resonance is not particularly broad, q ~ 6, the 
growth parameter is very large, \x ~ 0.3, and there is very efficient preheating in just a few 
oscillations, before backreaction sets in. At Hie same time, production of 4> and x particles is 
possible, although not as efficiently as for <rfc3 Note that preheating in this model could occur 
even in the absence of extra fields x-> simply due to the self-coupling of the symmetry-breaking 
field. 

Finally, in the case g 2 ~ A, typical of certain models of hybrid inflation in supergravity^ 
the two homogeneous fields oscillate with similar amplitudes and frequencies around the global 
minimum ((f) = 0, a = o"o). Since they are coupled, their frequencies and amplitudes vary rather 
chaotically and it takes many oscillations for their behaviour to stabilize around a periodic 
oscillation. During the chaotic motion, there can be no parametric resonance since this effect 
requires a periodic behaviour. By the time the oscillations become periodic, the amplitude has 
decreased so much that not even the narrow resonancecan be excited. Thus, for fh ~ H, there 
is no particle production in either of the three fieldsEj However, for fh 3> H, there are many 
oscillations in one Hubble time and eventually the motion becomes periodic while the amplitude 
of oscillations is still large. This results in a mild production of </> and a particles and an explosive 
production of x particles, if coupled only to the <f> field£3 



It should be noted that, in hybrid preheating, the equation that describes the instability 
growth of modes Xk is not exactly of the Mathieu type. Even if the x held couples only to one 
of the fields, <p or a, the fact that the two homogeneous fields follow coupled equations implies 
that the frequencies of oscillations of each field will vary strongly with time, unless the other 
field is fixed at its minimum. This does no preclude parametric resonance, it simply modifies 
the stability/instability chart. In some of the cases we studied^ the analysis with the Mathieu 
equation is a reasonably good approximation; in others, one has to compute the corresponding 
spectrum of instability bands. This is particularly important for the case g 2 ~ A, where both 
fields oscillate simultaneously in a rather chaotic way and their frequencies of oscillation depend 
explicitly on the value of the other field, which changes within one oscillation. 



3 Phenomenological consequences 

The processes by which the classical waves with large occupation numbers produced at pre- 
heating decohere and become an ensemble of relativistic particles in thermal equilibrium is still 
uncertain. Numerical lattice simulations have been performed to try to understand the pro- 
cess of backreaction and rescattering of those waves among themselves and with the inflaton 
backgrounds^ however, no one has yet been able to compute the final reheating temperature of a 
model in which preheating was important. It could be that, in the end, the reheating temperature 
is precisely the one computed from the usual perturbative analysis, T r h ~ 0.1\/TMp. However, 
the reheating temperature is not the only observable one can consider; in fact, preheating has 
opened the door to very exotic phenomena that might have novel experimental signatures with 
which to test our models of inflation. The interesting period is that after backreaction and before 
thermalization, where large numbers of particles (or waves) rescatter off themselves in a state 
far from equilibrium. Among these new phenomena, there are non-thermal phase transitions H 
with production of topological defects; generation of gravitational waves^ or production of heavy 
particles in a state far from equilibrium, which may help GUT baryogenesis Q or constitute the 
dark matter in our Universe! 

I will concentrate here on the production of a stochastic background of gravitational waves, 
which has specific signatures in the case of hybrid preheating. The idea is the following: the 
collisions among coherent waves of particles during rescattering radiates a reasonable fraction of 
energy, typicaly ~ 10~ 5 , in the form of gravitational waves. Their spectrum today depends on the 
details of rescattering and can be computed only numerically in lattice simulations! However, 
the low-frequency end of the spectrum, for wavelengths of order the size of the horizon at 
rescattering, can be computed analytically since in that case it is dominated by the gravitational 
bremsstrahlung associated with the scattering of x particles off the inflaton condensate, with 
the corresponding "evaporation" of inflaton particles. Taking into account that the occupation 
numbers of x particles are typically of order n,k(x) ~ 10 2 /iJ~ 2 at the end of rescatteringP and 
assuming that decoherence and thermalization occurred immediately after this stage, one can 
estimate the fraction of energy density in gravitational waves todays 

OgwMfc ~<W* > ( 5 ) 

where fi ra d^ 2 — 4 x ICU 5 and fig^u) = (p^ 1 dp gw / din u>)o. On the other hand, the wavelengths 
2n/k TS of the gravitational wave spectrum are redshifted to A ~ 0.5(M P i7 rs ) 1 / 2 A;- 1 todaytl In 
particular, the largest wavelength today corresponds to the size of the horizon at rescattering, or 
k TS ~ H rs , and therefore the minimum frequency in the spectrum is / m j n ~ 2 x 10 2 (Hj-s/TeV) 1 ^ 2 
Hz, which is right in the centre, for H IS ~ 1 TeV, of the range of frequencies detectable by LIGO. 
10 Hz < / < I0 4 Hz. In chaotic inflation, the fraction in gravitational waves today is typicallyta 
of order f2 gw /i 2 ~ 10~ 12 at the minimum frequency f m \ n ~ 10 6 Hz, which lies outside the range of 



frequencies and expected sensitivity of LIGO. The large rate of expansion at rescattering turns 
out to be a disadvantage for detecting gravitational waves from preheating in chaotic inflation. 
One would require a much lower rate, of order a few TeV, to be within the range of LIGO. 

Let us consider now a concrete hybrid inflation model with a low rate of expansion at the 
end of inflation, H ~ 10 2 TeV. Such a model could be consistent with COBE observations 
with very natural parameter values, g 2 ~ 0.01, A ~ 1,M ~ 10 12 GeV and m ~ 1 TeV. This 
model is of the first type, g 2 <C A, and has a large frequency of oscillations of the inflaton field, 
m ~ 10 e H. There will be explosive production of x particles even in the narrow resonance, 
q ~ 10 4 h 2 ~ 10 -3 , with a very small growth parameter /i ~ q/2. Backreaction in this model 
occurs at fht\ ~ (1/4^) ln(10 6 m(mti) 3 //if Mp) 1 / 2 ~ 33/q, which is still below one Hubble time, 
t\ ~ H^ 1 /SO, after the end of inflation. The occupation numbers of \ particles at this stage can 
be estimated^ as n/% (x) ~ 3 x 10 2 /if 2 q- 1 / 4 ~ 30 h 1 5/2 . Therefore, the estimate of Eq.(|) should 
be modified with an extra factor 10/ii in the denominator. Taking H TS ~ 10 TeV at the end of 
rescattering, the fraction of energy density in gravitational waves today is f2 gw /i 2 ~ 4 x 10~ 10 at 
the minimum frequency / m j n ~ 600 Hz, just in the appropriate range for detection by LIGO. 

We conclude that preheating in hybrid inflation has many features that differs from the 
usual chaotic or new inflation type of preheating, and makes it phenomenologically attractive: 
inflation may occur at a much lower scale, and still give the desired amplitude of temperature 
anisotropies in the CMB; there could be many oscillations of the inflaton field per Hubble time 
after inflation, which allows for a long-lived narrow resonance and very efficient preheating; and 
it is possible, in certain models, to generate a stochastic background of gravitational waves in a 
range accesible to observations. 



Acknowledgements 

I thank Andrei Linde for generous discussions and insights into the problems of reheating. 



References 



1. A.D. Linde, Particle Physics and Inflationary Cosmology (Harwood, Chur, 1990). 

2. L. Kofman, A.D. Linde and AA. Starobinsky, Phys. Rev. Lett. 73, 3159 (1994). 

3. L. Kofman, A.D. Linde and AA. Starobinsky, Phys. Rev. Lett. 76, 1011 (1996); 
I.I. Tkachev, Phys. Lett. B 376, 35 (1996). 

4. S.Yu. Khlebnikov and I.I. Tkachev, Phys. Rev. D 56, 653 (1997). 

5. E.W. Kolb, A.D. Linde and A. Riotto, Phys. Rev. Lett. 77, 4290 (1996). 

6. D.J.H. Chung, E.W. Kolb and A. Riotto, frep-ph/9802238] (1998). 

7. AA. Grib, S.G. Mamayev and V.M. Mostepanenko, Vacuum Quantum Effects in Strong 
Fields (Friedmann Laboratory Publishing, St. Petersburg, 1994). 

8. L. Kofman, |astro-pn/9605155| , |astro-ph/980222l| , |astro-ph/9802285| (1998). 

9. L. Kofman, A.D. Linde and AA. Starobinsky, Phys. Rev. D 56, 3258 (1997); 

P.B. Green, L. Kofman, A.D. Linde and AA. Starobinsky, Phys. Rev. D 56, 6175 (1997). 

10. J. Garcfa-Bellido and A.D. Linde, [hep-ph/9711360| , to appear in Phys. Rev. D (1998). 

11. A.D. Linde, Phys. Lett. B 259, 38 (1991); Phys. Rev. D 49, 748 (1994). 

12. C.L. Bennett et al, Astrophys. J. 464, LI (1996). 

13. L. Randall, M. Soljacic and A.H. Guth, Nucl. Phys. B 472, 377 (1996). 

14. J. Garcfa-Bellido, A.D. Linde and D. Wands, Phys. Rev. D 54, 6040 (1996). 

15. A.D. Linde and A. Riotto, Phys. Rev. D 56, 1841 (1997). 

16. S.Yu. Khlebnikov and I.I. Tkachev, Phys. Rev. Lett. 77, 219 (1996); 79, 1607 (1997); 
T. Prokopec and T.G. Roos, Phys. Rev. D 55, 3768 (1997); 

B.R. Greene, T. Prokopec and T.G. Roos, Phys. Rev. D 56, 6484 (1997). 



